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Motivated by recent experiments on electric transport through single molecules and quantum dots, 
we investigate a model for transport that allows for significant couphng between the electrons and a 
boson mode isolated on the molecule or dot. We focus our attention on the temperature dependent 
properties of the transport. In the Holstein picture for polaronic transport in molecular crystals the 
temperature dependence of the conductivity exhibits a crossover from coherent (band) to incoherent 
(hopping) transport. Here, the temperature dependence of the differential conductance on resonance 
does not show such a crossover, but is mostly determined by the lifetime of the resonant level on 
the molecule or dot. 

PACS numbers: 73.23.-b,73.63.-b,71.38.Fp 



I. INTRODUCTION 

In recent years there has been a growing mierest in 
electrical transport through single molep^lesErcfuiJ and 
single electronic levels in quantum dotsijO'LJ Some molec- 
ular devices exhibit switching behavior with large on-ofF 
ratiosu increasing thej-|motivation to construct molecu- 
lar electronic devices.Q In some cases it has been found 
that the transport is Wiite temperature dependentEl and 
it has been suggestedcl that this is due to the presence 
of low energy boson modes, such as internal rotations, 
which couple strongly to the molecular electronic. iStatcs, 
and can easily be excited by small temperaturesJj'EEl In a 
similar vein, in double quantum dots it has been found 
that there arajacoustic phonons which couple strongly to 
the electronsElEl. 

Some experimental values for the phonon energy have 
been estimated in various papers. In Table ^ we give some 
numbers for reference. We see that the boson (usually 



TABLE I: Typical values for parameters taken from experi- 
ment, huio is the boson energy and F is the line width due 
to coupling to the leads (defined below) of the resonant level 
on the molecule or dot. Imax is the maximal current driven 
through the system. 

System ^ F Imax 



2 quantum dots.i— ■ 40 fieV 0.2 ^eV 3 pA 
2 quantxHjpi dotajtil 30 peV 1 /leV 5 pA 
moleculetltl 



Ceo moleculen 



3 meV 
5 meV 



1 nA 
0.1 nA 



phonon) frequency in these systems is quite small, corre- 
sponding to temperatures in the range 0.|5-|- 50 Kelvin. 
In addition there was a recent proposaM to consider 
transport through a quantum dot to a carbon nanotube 
cantilever with a resonant frequency of the order of 100 
MHz, corresponding to a phonon energy of 0.4 /icV. If 
the electron-phonon coupling is sufficiently large pola- 
ronic transport might be important for these systems. 
When the electron tunnels through it can absorb or emit 
bosons, thus altering its energy and the current. If the 



temperature is much larger than the boson energy, there 
are many bosons available for absorption and this might 
heavily influence the current. 

In 1959 HolsteinEi predicted that for a periodic 
one-dimensional molecular crystal with strong electron- 
phonon coupling there should be a crossover from co- 
herent (band) to incoherent (hopping) transport with 
increasing temperature. When increasing the temper- 
ature the effective bandwidth becomes narrower, this 
gives rise to a decrease in coherent transport. In con- 
trast, increasing temperature means that more and more 
phonons are activated and we are in a regime where 
phonon assisted inter-site tunneling starts to contribute 
to the conductivity. This coherent-incoherent crossover 
is believed to have been observed for the first time quite 
recently in single crystals of pentacene.Ej One aim of 
this paper is to see whether a similar crossover should 
be seen in polaronic transport through molecules and 
quantum dots. This might be expected because of the 
mathematical similarity between the models for periodic 
systems and the resonant tunneling case. We might 
expect the tunneling amplitude between the leads and 
dot to be reduced by polaronic effects, thereby reduc- 
ing the coherent part of the conductivity. When increas- 
ing the temperature the electrons can tunnel with bo- 
son assisted transport that enhances the tunneling, pos- 
sibly leading to a crossover behavior. There have been 
many theoretical investigations of the. ff[f<]t' pf phonons 
on the tr ansupr |Ut|hmiigh moleculeaaa'OE^nlZl and quan- 
tum dotMBBme, but none of them focuses on the 

temperature dependence of the current. The purpose of 
this paper is to clarify this aspect of the transport. Li 
et. alc3 included a Hubbard term, but did not consider 
multi-phonon contributions. In a recent paper Emberly 
and KirczenowllZI made a thorough analysis of conduc- 
tance through a molecular wire. A set of self-consistent 
equations where set up and solved to give the distribution 
functions in the leads and molecule, and then transmis- 
sion probabilities were calculated. However, the temper- 
ature dependence is not addressed in that paper. 

In this paper we perform the analysis for the simplest 
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possible case, where the electrons interact with a single 
optical boson localized on the dot or molecule. We antic- 
ipate that this is sufficient to illustrate the main physics 
in the more complicated case of many bosons, such as 
acoustic phonons. In order to obtain analytical results 
we have to assume that the coupling to the leads is small 
and the energy level in the dot or nmLeciile is not too 
close to the Fermi energy in the leadsO'E£l By assuming 
that the coupling to the leads is small we can calculate 
the effects from the bosons locally on the molecule/dot 
and then assume that the effect on the leads from the 
bosons is negligible. This enables us to use well-known 
results from mesoscopic transport theory. The bosons 
are possibly most often phonons, but since the theory 
will look identical (assuming linear couping) for different 
types of bosons (phonons, magnons, charge oscillations) 
we will simply refer to "bosons". Even in photon as- 
sisted tunneling through quantum dots side bands have 
been observed when tuning the photon energy.EEl In sec- 
tion II we will define the model we use and in section III, 
we discuss the approximations we have to make. Differ- 
ent limits for the current are derived in section IV, and 
in section V we discuss the differential conductivity. 



II. CURRENT THROUGH A LEVEL COUPLED 
TO A LOCAL BOSON MODE 



is given by 



2-3, 



(1) 



where 



ki 



'Hl-2 = ^ ^(CkjC2 + h.c), 
ki 



Here eo is the energy of the level in the dot/molecule 
and t is the energy associated with hopping onto/off the 
dot. The electronic dispersion in the leads are given by 
and • M is the coupling to the local boson mode 
with energy hwQ. We disregard the spin dependence for 
simplicity. 

First we make a unitary transformation to diagonalize 
the Hamiltonian 7^2 . The price we pay for this is that ex- 
tra operators attach to the tunneling term in the Hamil- 
tonian. The transformation is 7^2 = e^Ti.2e^^, where 
S = c1'c^(a^ - a). This gives us 



7^2 



huoCL'a — /S.a2C2, 



(2) 



We consider the simplest possible model Hamiltonian 
and neglect the spin degree of freedom and any effects 
of electron-electron interactions. The system we study 
consists of the individual entities (left lead, molecule or 
quantum dot, and right lead) coupled via tunneling. We 
assume that we are dealing with a resonant tunneling sit- 
uation, but the states in the dot (or molecule or any single 
level system) couples to some boson mode with character- 
istic frequency wq, as shown in Fig. |l|. The Hamiltonian 



where 



'Hl-2 7^2-3 



eV/2 




eV/2 



FIG. 1: Tunneling through a system with one level. The 
dashed lines indicate the bosonic satellites (see text). The 
Fermi energy in the leads is chosen to be zero. The electrons 
has to tunnel through the barriers, and can absorb or emit 
bosons in the process, corresponding to the lines below and 
above the central resonance respectively. The Hamiltonian 
given in Eq.(|^) contains terms describing the different parts 
of the system. 



A = - eo. 



(3) 



When the central system is a quantum dot eo (and thus 
A) can be adjusted by applying a gate voltage. After 
the transformation the tunneling part of the Hamiltonian 
becomes 



7^i-2 = X!*(4i^2^ + ^-^-)' 

ki 

H2-3^Y.*{cl^''^X + h.C.), 



where 



X = exp 



— {a -a?) 



(4) 



(5) 



The X-factors can be absorbed into a 
electron creation/annihilation operator 



renormalized 
in region 2, 

so that we are left with the usual resonant tun- 
neling Hamiltonian except that the Greens function 
for the electrons on the molecule/dot has an addi- 
tional complication. (Trc(r)ct(O)) (r^c(T)ct(O)) = 
{Trc{T)c'^{0)){TrX{T)X'^{Q)). A formula for the curxepl 
can be derived using a Landauer-Biittiker approach. EZIEEl 
First we calculate the current from the left lead onto 
the dot from the rate of change of particles in the left 
lead. A similar expression for the current from the dot 
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to the right lead is derived and the total current through 
the system is obtained by combining these two formulas. 
The derivation is presented in detail in Ref. 27 and |2|. 
The result is that the current is given by 



I{V) 



2e 
T 



*[/i(e)-/3(e)]Im[tr(rG2(6))]. (6) 



The applied voltage across the system is V and it enters 
the two Fermi functions (the equilibrium Fermi level of 
the leads is chosen to be zero) /i(e) — /(e — eV/2) and 
/3(e) = /(e + ey/2). Further, 6*2(6) is the Green function 
for the quantum dot including all effects from the boson 
system and the tunneling to the leads. The parameter F 
is 



TlFg 



(7) 



where ri(3) = 27rt^_Di(3)(e), -Di(3) is the density of states 
(DOS) in the left (right) lead. Fi(3) is the width of the 
central resonance due to the tunneling to the left (Fi) and 
right (r3) lead. The total width of the local resonance, 
is the sum of the two, r2 = Fi + F3. 
For convenience we introduce the dimensionless param- 
eters, 
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M ■ 

hujo 

F 



We emphasize that there are many different energy scales 
associated with the system: fc^T, eV, tvjjQ, F, M, and eq- 
The relative sizes of these energy scales have a significant 
effect on the current through the system and what ap- 
proximations can be made in evaluating it. 

The electrons will deposit/absorb energy from the 
bosonic system that has to be carried away/supplied. 
Therefore a question arises about how to define the tem- 
perature, particularly of the molecule or dot. We assume 
that the molecule/dot is in equilibrium with a bath and 
that the tunneling rate is small so that the system relaxes 
to the initial state after each tunneling event. In a quan- 
tum dot the bath can be the substrate that the quantum 
dot is manufactured on. For a molecule a surrounding 
cooling liquidEJ can play the role of the bath. Otherwise, 
we have to assume that the deposited or absorbed energy 
is transferred to/from the molecule via the leads. As far 
as we are aware, this assumption is also (implicitly) made 
in all other theoretical work on this subject. 



III. APPROXIMATE EVALUATION OF THE 
GREENS FUNCTION G2(e) 

To be able to use Eq.(^) we have to calculate the lo- 
cal Green function, 6*2 (e)- Due to the coupling to the 
leads finding G2(e) is a highly non-trivial problem in 



many-body theoryBd-EilEa It is comparable in difficulty 
to the Kondo problem because of the possibility of non- 
perturbative effects. This is true even in equilibrium (i.e., 
in the absence of a bias, V = Q). A recent study was 
made of a similar Hamiltonian_Iwith spin) using the nu- 
merical renormalization group. lj We are interested in the 
non-equilibrium case where there is a bias. In order to 
simplify the analysis we have to rely on approximations, 
and the result will depend on how the X-operators from 
Eq.(||) are decoupled. One alternative is to assume that 
the coupling to the leads is small, F1-I-F3 ^ A, this is the 
approach taken here. This approximation is justified for 
small currents, as is the case in the systems considered 
here. If we were to include the effect on the leads from 
the bosons on the molecule/dot there would be a nar- 
rowing effect on F. HewsotLand Newns used variational 
and perturbation methodstj to show that this narrowing 
only takes place if the following conditions apply: 

— > F, gi > 1, hujQ > Fe^fi 
nujQ 

hujo > |A|. 

The conditions on the first line means that the electron- 
boson coupling has to be large enough to form a polaron. 
The last requirement on the first line means that individ- 
ual boson satellites can be distinguished from each other. 
The second line tells us that if the level and boson satel- 
lites are too far from the Fermi level it is energetically 
unlikely to have virtual boson excitations, thus the leads 
are unaffected by the bosons. The narrowing is approxi- 
mately given by 



t 



te 



-Si(l/2+ni5) 



where is the Bose function. 



1 



(8) 



(9) 



and f3 = l/ksT. 

The above considerations apply to equilibrium iV — 0) 
whereas we are interested in the non-equilibrium situa- 
tion of a finite bias, and particularly the resonant tun- 
neling case where one of the leads' Fermi level is close to 
the dot/molecule level {eV = ±A/2). In that case the 
narrowing of the level width due to that lead (but not 
due to the second lead) may occur, e.g.. 



r2 = Fi+F3^Fi+F3e 



-gi(l+2nB) 



(10) 



If Fi ~ F3 this will lead to some quantitative but no sig- 
nificant qualitative changes in the current-voltage char- 
acteristics and so we will not consider them further. 

We treat the leads as unaffected by the bosons, i.e., no 
narrowing of the bands in the leads. This means that we 
ignore the averages of the X-operators that appear in the 
tunneling part of the Hamiltonian, Eq.(^, the justifica- 
tion for this is given above. Below we will also assume 
that the leads give rise to a flat, energy independent. 
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density of states. This is sometimes called the wide band 
limit .E3 Otherwise T would be energy dependent. The 
quantum dot Green function calculated using these ap- 
proximations is 



G2(t) = -^e(^)e(^^-^=/2)t/Bg-^•(t). 



(11) 



The factor e ^ii' is due to the coupling to the boson and 
can be writtenEa 



g-*(t) ^ g-ffl(l + 2«B) 



, (12) 



where /; denotes a modified Bessel function. 

We Fourier transform the Green function and get an 
expression for the total current 



~h 

CO 



^-lhuJo/3/2 



de [/i(e)-/3(e)]e-^^<'+'""^ 

) 

Ii \^2gi^/nB{l + nB) 

L — — OQ 

T1+T3 

(e + A + /ftwo)2 + ^^^^^ ■ 



(13) 



We can interpret Im [tr(rG2)], in Eq.(H), as the trans- 
mission coefficient for the tunneling. We plot this in 
Fig. ^ for a certain choice of parameters. The resonances 
to the left of (e -I- A)/huJo corresponds to absorption of 
bosons, and the ones to the right to emission of bosons. 
The middle line can be identified as the so called zero- 
boson transition. The width of each satellite depends on 
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FIG. 2: Transmission coefficient, Im[tr(rG2)], as a function 
of the energy, for three different temperatures. The satellites 
are due to the boson modes. g\ = {M/hojoY = 0.5 and 
(/2 = (r/fttJo)^ = 0.09. The vertical axis is normalized to the 
highest peak in the plot. 

r2 directly. When increasing the temperature the satel- 
lites increase in amplitude, indicating that it is easier to 
emit /absorb bosons. The asymmetry between negative 



and positive energies is due to the factor Q-'-f^oP/'i_ This 
is a due to the fact that at low temperatures there are 
no available bosons to absorb. 

In Fig. ^ we plot the current as a function of volt- 
age using Eq.dl3) for a set of parameters. It show steps 




FIG. 3: Current as a function of the applied voltage for dif- 
ferent choice of couphng strengths. We set ksT — 0.1ftti;o and 
A = hu}o/2. 



indicating that more and more satellites participate in 
conducting electrons. Note that the steps in Fig. || occur 
every second Hluq. This is simply because the satellites 
are positioned equidistant on each side of the central res- 
onance, we have to increase the voltage by 2hujo in order 
to cover the satellite. The first satellite starts to con- 
tribute to the current when eV = 2 A. A decrease of r2 
{g2 decrease) results in sharper steps, and a decrease in 
the amplitude of the current. When F :s> huiQ (large 32) 
the step structure disappears. Increasing the tempera- 
ture results in the step structure being washed out to a 
smooth curve. 

When increasing M, the amplitude of the current drops 

due a decrease of the factor e i'^^o) Eq.(p^. 
Increasing the temperature has the same effect. With- 
out any coupling to the boson {M, gi—0) we get a single 
resonant level without any satellites. This can be seen in 
Fig. 1^ where we plot the current as a function of eo, the 
location of the energy level in the dot or molecule. The 
application of a gate voltage in a quantum dot -would 
be equivalent to changing the level eo (or A) .00 We 
see a shoulder developing corresponding to the first bo- 
son satellite. A similar effect has been seen in a double 
quantum dot system.cl The absence of a boson absorption 
peak in Fig. ^ is due to the low temperature, this comes 
from the factor g"'""''/^. If we increased the tempera- 
ture, or the electron-boson coupling, enough there would 
be more side bands visible. 
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FIG. 4: Current as a function of the location of the energy 
level in a quantum dot when bosons are present {gi — 0.1) 
and absent {gi = 0). ksT = O.OShtoo and we put g2 = 0.5. 
eV is set to 0.2huJo so that we only scan a small region around 
eo. Parameters are taken from ref. ^. We only see the boson 
emission satellite due to the low temperature and the small 
electron-boson coupling. 



2. T = 

If the temperature goes to zero we can approximate the 
Fermi functions with step functions. Then, the integral 
over e can be performed and the result is 



lim I(V) 



2er 

IT 



tan 



1 f eV-2eo 



tan 



1 f eV + 2eo 
(16) 



Further, if eV and 2eo is small compared to Fi + Fa we 
can use the property that tan^^(a;) ~ x, and we get 



lim 

(ey,2eo)Ari+r3)<i 



4e^F 



/^(ri+Fa) 



V, 



(17) 



i.e., a linear regime at low voltages. If, on the other hand, 
we take 1/ ^ oo in Eq.(p^) we get 



lim / = — . 

V^oo 



(18) 



IV. LIMITING BEHAVIOR FOR THE 
CURRENT 



This means that the whole resonant level contributes 
maximal to the current. 



In order to better understand the influence from the 
bosons on the current. Let us now have a look at the 
current in some limits. 



B. A/ / 

1. eV">fcsr, ficjo 



A. M = 

If we put the coupling between the boson and the elec- 
trons to zero, we get: 



TlFg 



(14) 



This would correspond to resonant tunneling without any 
bosons. 



1. kBT:$>T,eV 



In this case we get the same limit as in Eq.(|^) even 
if AI 7^ from Eq.(|l^). This can be seen in Fig.^ where 
all curves tend to the same value at large V. If we have 
that eV ^ AjksT we can replace /i(e) — /3(e) by a 
factor 1, and the integral would extend between —eV/2 
and eV/2. But since eV is greater than all other energies 
we extend the integral from —00 to 00. The integral gives 
a contribution tt. All parts coming from the boson gives 
1 and we again have the limit 



/(el/ > A, /cbT) 



eF 



(19) 



This limit can be seen in Fig. ^ where all curves tend to 
the same limit at high applied voltage. 



In this limit Eq. (|4|) reduces to the linear response ex- 
pression 



ttF 



!_ _ 4e2 

y~ h kBT cosh' 



(15) 



A similar form was used by Qin et. aB to fit their exper- 
imental data. 



2. ksT < hujo 

Let us now investigate the limit ksT <^ Hloq. In this 
limit (corresponding to low temperatures) we can ap- 



proximate the Bose function as 



< 1. 



All terms corresponding to positive I vanishes. This is a 
result of the physical fact that positive I corresponds to 
boson absorption but at T = there are no bosons. The 
Bessel function can be approximated as Ii{z) ^ jf{z/2)' 



6 



when z — > 0. Then we get that the current becomes 



T 2er 

ikBT<^hwa = j^e 

lil 



de[fi{e)-Me)] 



Ti+r., 



1= 



(e + A + Ihujo^ + 



(20) 



3. ksT > hujo 

For high temperatures we approximate the Bose func- 
tion as — kBT/hujQ. The argument in the Bessel 
function is large and we can use the property 



z > 1. 



(21) 



Using this, the current becomes 

3^91 4fer,T 



2er 



1— — OC 



{e + A + lhwoY 



4 



de[h{e)- 

(22) 




FIG. 5: Failure of the saddle point approximation. Current 
calculated in two different ways: the full lines were obtained 
using the exact result (Eq.(y_3|)) and the dashed lines using 
the saddle point approximation (Eq.(^5|)) in the expression 
for the current. We see that the saddle approximation does 
not reproduce the full expression for the current. Here we set 
ksT = G.lhwo and A = hujo/2. 



DIFFERENTIAL CONDUCTANCE 



C. Saddle point approximation 

If 52 ^ 1 (i-e., r ^ ?iwo) we can evaluate the current 
using a saddle-point approximation similar to that used 
previously in Ref. ^ and [l^. The exponential factor, 
e-*(*) EE (X(t)Xt(O)), in Eq^) can be written as 



-9l[(nB + l)(l-e"""«')+«f3(l-e""0*)] 



(23) 



We approximate the exponential function in the expo- 
nent, ^ 1 + z + and we get: 



G2{t) 



pi{A+gihuJo)t/h-r2t/2h~gi/2{l+2nB)iuiQty 



(24) 



Let us assume that we can neglect the term linear in t in 
the exponent compared to the quadratic one, i.e., gi 3> 
g2- We Fourier transform the resulting Green function 
and get that the relevant factor entering Eq. (ph becomes 



Ini[G2(e)] 



exp 



' 2gi{hLUoy^l+2nB} 



ujo\/gi{r+2nB) 



(25) 



This approximation gives a broad Gaussian line shape 
covering all the boson satellites. This is in contrast to 
the individual boson satellites shown in Fig. ||. Using the 
saddle point approximation would give a Gaussian line 
shape in /(eo)i whereas a Lorentzian line shape occurs in 
the regime, ksT ^ F ^ Swq, illustrated in Fig. |^. 

The current using Eq.(p5|) is plotted in Fig.|^ and com- 
pared to the full expression, Eq. (p^) . In this figure we 
can clearly see that the saddle point approximation can- 
not reproduce the actual current. Only for a small range 
of bias voltages, for low temperature and large coupling 
is there an agreement. 



The differential conductance, defined by 

dl 



(26) 



more clearly reveals the effect of the bosons. In general 
this is given by 



de{h{e) [l-A(e)] + /3(e) [l-/3(e)]} 

) 

OO 

3l(l+2n.) ^ [251^^5(1+^5) 

/=- 

lhiAjol3/2 



2 



{e + A + Ihwo)'^ + 



(27) 



Later we will set eV = 2A, which corresponds to reso- 
nant transport through the zero phonon feature. If we 
let the temperature go to zero in this expression we can 
approximate the Fermi functions together with the tem- 
perature as a delta function, /3n/(e)[l — ?^/(e)] ^ 
and again only negative / contributes, corresponding to 
emission of bosons, and we get 



[C{V)]t^o = 



e^FiFg 
2h ' 



W 



1 



(A + Ih^o + eU/2)2 + ^^1^ 
1 



(A + Ihiuo - eU/2)2 + (£i+E3)! 



(28) 
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We define 



{a. 



[C{eV = 2A)]t- 



(29) 



For the particular case A 3> fiujo, gif^o, this simphfies 
to (Cres)" — ^ Ti+Fa ' showing how polaronic eff'ects 
reduce the diff'erential conductance. 

In Fig. II we plot the differential conductance as a func- 
tion of the applied voltage for difl'erent values of temper- 
ature and coupling parameters. The peak at eV = 2A 




FIG. 6: Differential conductance as a function of applied volt- 
age when changing the temperature (upper left), electron- 
boson couphng g\ (upper right and lower right) and the level 
width g2 (lower left). A = fttJo/S. At eV = 2A it has a 
maxima for moderate couplings gi < 1. To obtain a maximal 
signal it would be desirable to perform the experiments at 
this value. 

correspond to the zero-boson peak, and in the consecu- 
tive peaks one/two/three. . . , and so on, bosons are emit- 
ted or absorbed. As seen in this figure increasing the 
temperature, or the level widths, drastically affects the 
shape of the differential conductance. 

In Fig. we plot the differential conductance on reso- 
nance with the zero phonon line as a function of temper- 
ature for a range of parameters. In this figure we see that 
the differential conductance generally decreases with in- 
creasing temperature, in contrast to the non-monotonic 
dependence found by Holsteinl3 for periodic molecular 
crystals. The corresponding crossover behavior does not 
occur for transport through molecules/quantum dots, 
since this would be indicated by an upturn in Fig. ^ when 
increasing the temperature. The absence of a crossover 
can also be seen by looking at Fig. ^ from that the slope 
at eV = 2 A (i.e., the differential conductance) is almost 
constant when changing gi from to 0.5. If we were to 
calculate the differential conductance in the limit when 
ksT 3> fiwo from Eq.(p^ wc see that the temperature 
dependence of the differential conductance is governed by 



the pre-factor ■ 



and this is a strictly decreas- 



ing function of the temperature, for reasonably values of 
gi. Thus, there will never be an upturn in the differential 



a 0.6 



U 



U 







1 






8 


=0.0 


g,=0.02 






-- g 


= 1.0 


g,=0.02 






.... g 


=0.0 


g,=l-0 






■-■ g 


=2.0 


g,=l-0 




















1,1,1,1,1 




FIG. 7; Differential conductance at the resonance as a func- 
tion of temperature for different parameters. The upper graph 
shows that, for moderate couplings gi < 1, the differential 
conductance is ahnost unaltered by the presence of the bosons. 
The lower graph shows that the differential conductance, for 
moderate couplings f/i < 1, is determined by the parameter 
ksT /T. The plots were made assuming a constant DOS in 
the leads, A = and eV = 0. 



conductance when increasing the temperature. This gen- 
eral behavior is not changed when gi is changed. Even 
an increased applied voltage, meaning that more boson 
satellites contributes, was not able to induce a crossover. 
However, changing F does alter the amplitude of the dif- 
ferential conductance, as seen in Fig. |^. 

As mentioned above the temperature behavior is dom- 
inated by F. If we put M = in Eq.(|27|) we can write 
the differential conductance as 



res\M-*0 



h ksT J_ 



dy[f'{y) + f'{y 



2eo 



1 



y2 + {T/kBT)^ 
(30) 



where f = (Fi -t- F3)/2. If we now take eo = or 
ksT Co we will have that the differential conductance 
is a universal function of F/fc^T, i.e. 



[C, 



reslM^O 



F{f/kBT). 



(31) 



This can be seen in the lower graph in Fig. |^, where the 
two graphs for (?i = (but g2 = 0.02 and 2.0 respectively) 
collapse on the same line when the temperature axis is 
rescaled. 



VI. CONCLUSION 

In conclusion we see that the polaronic transport 
through a single molecule or quantum dot does not 
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clearly exhibit the crossover from coherent to incoherent 
transport expected for the case, of a periodic molecular 
crystal considered by Holsteinlla. The general behavior 
of the temperature dependence of the differential con- 
ductance is in large unaffected by the presence of the 
bosons. The temperature dependence is mostly deter- 
mined by the linewidth (due to coupling to the leads) of 
the resonant energy level. The bosons produce side bands 
corresponding to absorption and emission of bosons. We 



also stressed that because of the interaction of the po- 
laron on the dot or molecule with the leads there are po- 
tentiaiy .same very interesting problems in many-body 
physicsEjCj to be explored in the model system we have 
considered. 
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